A new polymorphic material? Structural degeneracy of ZrMn2 
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Based on density functional calculations, we propose that ZrMn2 is a polymorphic material. We predict that 
at low temperatures the cubic C15, and the hexagonal C14 and C36 structures of the Laves phase compound 
ZrMn2 are nearly equally stable within 0.3 kJmoP^ or 30 K. This degeneracy occurs when the Mn atoms mag- 
netize spontaneously in a ferromagnetic arrangement forming the states of lowest energy. From the temperature 
dependent free energies at T ~ 160K we predict a transition from the most stable C15 to the C14 structure, 
which is the experimentally observed structure at elevated temperatures. 

PACS numbers: 61.50.Lt, 71.15.Nc, 61.66.Fn, 75.50.-y 
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Polymorphism of a solid material is its property to exist in 
several or even many crystal modifications which are called 
poly types. A very famous example of a polymorphic mate- 
rial is SiC |1] for which the polytypes are distinguished by 
the variations of stacking of SiC bylayers. ZnS, as well, ap- 
pears in many polytypes |2]. First principles calculations (as 
usual referring to T=OK) on SiC 1 3] and on ZnS and related II- 
VI compounds |4] revealed that the occurence of polytypism 
must be due to the energetical degeneracy of structures at very 
low temperatures. For Laves phase compounds 1 5] so far only 
high temperature polymorphism is known from experiments 
on e.g. TiCr2, ZrCr2, and HfCr2 0]. The polytypes are the 
most common structures of Laves phases compounds, namely 
the densely packed hexagonal C14 and C36, and the cubic 
C15 crystal structures which are distinguished by the stacking 
of layers and blocks as sketched in Fig. ^for ZrMn2. Search- 
ing for the ground state structure of ZrMn2 we performed den- 
sity functional calculations for the C14, C15 and C36 struc- 
tures finding them energetically nearly degenerate at very low 
temperatures. This new finding, the first case of low tempera- 
ture polymorphism of a Laves phase compound, is the subject 
of our paper. According to our results magnetism plays an 
important role in the polymorphism of ZrMn2. Bulk samples 
of ZrMn2 were fabricated so far only at high temperatures of 
about 870 K at which the C14 structure was observed |6, 7]. 
Our calculations confirm the stabilization of this high temper- 
ature phase due to the free energy of the lattice vibrations. 

We predict that at least the C14, C15 and C36 structures 
are nearly equally stable at very low temperatures. All three 
phases are ferromagnetic and their structure dependent gain 
in magnetization energy is at just the right size to cause the 
degeneracy. The CI 5 structure, which is the most stable one 
at T=OK, is predicted to transform at Ttr ~ 160K to the C14 
structure as observed experimentally at elevated temperatures. 

Another type of structural degeneracy at low temperatures 
might occur due to martensitic instabilities by which ener- 
getically closely related structures are transformed into each 
other because of elastic or vibrational instabilities i3>@l- For 
ZrMn2, however, our calculations reveal no martensitic in- 



stability in the analysis of elastic and vibrational properties 
which strengthens our argument that ZrMn2 is a truly poly- 
morphic material. 

Fig. n schematically shows the construction of the C14, 
C36 and CI 5 crystal structures in terms of stacking of blocks 
consisting of Zr and Mn layers. As is visualized by the figure 
and described in the caption, the stacking is such that blocks 
can only be combined when the layers are arranged appropri- 
ately. Following this recipe one might construct longer ranged 
stacking with larger unit cells 01 which possibly might lead 
to a large number of polytypes. 

For the calculation of ground state properties of ZrMn2 
in the C14, C15 and C36 crystal structures we applied den- 
sity functional theory (DFT) by means of the plane wave Vi- 
enna Ab initio Simulation Package (VASP^ Ll 1.1 which -in its 
projector augmented wave formulation - is one of the most 
precise methods for calculating the energetics and electronic 
structure of solid matter within periodic boundary conditions. 
For the many-body exchange-correlation interaction the gen- 
eralized gradient approximation of Perdew and Wang 1 12] in 
combination with the approach of Ref. L13iil for spin polar- 
ization was chosen. The calculations are free from any em- 
pirical parameters. Lattice parameters as well as the atomic 
coordinates are determined by minimizing total energies and 
atomic forces. Care was taken that for each structure a suf- 
ficient number of k points for the Brillouin zone integration 
was chosen. Ferromagnetic as well as some selected antifer- 
romagnetic spin arrangements were considered. Details will 
be described elsewhere |14]. In order to take into account 
the temperature dependent effects of lattice vibrations, the 
phonon dispersion and density of states were calculated by 
the so-called direct approach 1 15] utilizing data of VASP cal- 
culations. Furthermore, elastic constants for the C14 and C15 
structure were determined from DFT energy densities corre- 
sponding to suitable distortions. The energy of formation was 
obtained from the difference of total energies for the com- 
pound and the pure phases fldll . 

Figure 2 contains the essential information of our study. 
Viewing the nonmagnetic (NM) results, clearly the energy 



2 



O Zr atom 
' Mn atom 

(a) C14 [0001] 



(b)C15 |ni| 

mil 





;Z20=:fib==0 



FIG. 1: Schematic stacking of Zr and Mn layers for the hexago- 
nal C14, C36 and cubic CI 5 crystal structures along the directions 
[0001] and [111], respectively. C14 (panel a) is build by the two 
blocks Bl and Bl'. Bl consists of layers Zl with one Zr atom per 
layer unit cell, and of two planes Ml(l) and M2(3). Ml(l) has 1 Mn 
atom per layer unit cell, whereas the unit cell of layer M2(3) con- 
tains three Mn atoms arranged in triangles typical for the structures 
of Laves phases. Primed objects (blocks or planes) are obtained when 
mirroring the corresponding unprimed object by the mirror plane 
shown in the insert (left top). The C15 structure (panel b) is built 
by the three blocks B2, Bl', and B3. Block B2 contains the Mn layer 
M3(l), and in block B3 the layer Z2 has Zr atoms of weight 1/4 at 
the corners. The layer M4(3) comprises 3 Mn atoms per layer cell 
(Mn atoms on boundary lines have weight 1/2, the atom in the center 
counts fully). The C36 structure (panel c) is constructed by stacking 
of blocks B2', B3, B2, andB3'. 



of the C14 structure is lowest for all volumes. At equilib- 
rium volume Vo corresponding to the minimum of the energy 
curves C14 is favoured by about 1.5 kJ mol^^ compared to 
C36, and by the substantially larger amount of 5.3 kJ mol^^ 
compared to C15. The important role of magnetism is al- 
ready illustrated by the results for the most stable cases of 
selected antiferromagnetic ( AF) spin arrangements 1 14^ .16il . 
For all structures, the minimum energies are lowered, lying 
now much closer together as compared to the NM case, with 
C36 as the favored structure. Also noticeable is that the equi- 
librium volumes for C14 and C15 are quite different express- 
ing different magneto-volume effects. The most striking result 
is presented by the ferromagnetic (FM) data. Again, the en- 
ergy is lowered for all structures when compared to AF, but 
now C14, C15, and C36 nearly coincide. By the very small 
energy difference of 0.3 kJ mol^^ (corresponding to « 30K) 
the cubic CI 5 phase is now more stable than the two hexag- 
onal cases. At larger volumes, differences are getting slightly 
larger, whereas at smaller volumes (i.e. under pressure) the 
energy curves of all three structures are becoming undistin- 
guishable. 

The interpretation of these results is that for the crystal- 
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FIG. 2: Density functional results for ZrMn2 . Left panel: Calculated 
energies Udft{V) as a function of volume V for the nonmagnetic 
case (NM, asterisks and full black lines), the most stable antiferro- 
magnetic (AF, triangles and dashed lines), and ferromagnetic (FM, 
dots and full lines) spin orderings for the C14, C15 and C36 struc- 
tures of ZrMn2. Volume per formula unit (3 atoms per unit cell). For 
each marked point of the curves all structural parameters are fully 
relaxed resulting in the minimum of total energy for a given fixed 
volume. Right panel: Energies Udft at equilibrium volume Vo at 
minimum of the total energies of all the cases in the left panel. The 
zero of energy corresponds to lowest calculated energy which was 
obtained for C15 FM. 



lization process at least at very low temperatures there is a 
choice between at least three structures distinguishable by the 
stacking sequences according to Fig. ^ If one wants to fab- 
ricate low temperature poly types of ZrMn2, molecular beam 
epitaxy seems to be the only choice. Standard procedures 
of bulk crystal growth are rather unsuitable for the creation 
of low temperature polytypes, because at high temperatures 
ZrMn2 crystallizes only in the C14 structure as observed at 
about 870K 0). On cooling down the C14 structure might be 
frozen in as a metastable state. The thermodynamical equilib- 
rium, which according to our prediction is the CI 5 structure, 
might be never reached in a finite time. In addition, no sim- 
ple deformation path can be imagined to deform C14 to C15 
or C36, making a structural transition even less probable. At 
present, no unambigous experimental structural information 
for low temperatures is available. 

ForZrMn2, Stoner's theory predicts the occurence of band 
magnetism already from the NM density of electronic states at 
Fermi energy. The energetical degeneracy of structures then 
occurs for the feiTomagnetic (FM) case. The most important 
indication why this happens is revealed by the right panel of 
Fig. 12] The FM energies for the coiTesponding Vq are al- 
most lying on one line quite in contrast to the NM case for 
which the C15 energy is rather high compared to C14. By 
spontaneous FM spin polarization all three structures gain dif- 
ferent magnetization energies defined by Udft{FM; Vq) — 



3 



Udft{NM; Vq), which are just of the right size to make the 
structures degenerate. The largest magnetization energy of 
9.2 kJ mol^^ is gained in the cubic C15 structure. The de- 
generacy clearly is due to an electronic effect, since the struc- 
tural changes, when relaxing geometries for all the magnetic 
states, have only insignificant influence on the energy differ- 
ences. The average magnetic moment of Mn in all cases is 
about 0.9 /i_B. It should be noted, that for the hexagonal cases, 
not all Mn atoms are equivalent, therefore also the magnetic 
moments differ significantly. The average magnetic moments, 
however, agree very well between the structures. For the C14 
structure the magnetic moments of the two non equivalent Mn 
atoms are 1.24 and 0.60 /is, for C15 the Mn moment is 0.90 
/is, and for the C36 case three non equivalent Mn atoms ap- 
pear with moments of 1.14, 0.92 and 0.71 /.ib- In contrast to 
other compounds with long range atomic stacking sequences 
flTil . which are non magnetic, for ZrMn2 the ferromagnetic 
spin ordering for the energetical degeneracy of suitable atomic 
stackings is crucial. 

Concerning antiferromagnetic spin orderings, it might be 
possible that more complicated AF structures (including non 
collinear spin arrangements) might lower the AF total energy. 
However, also due to the rather small magnetic moment we 
would not expect any substantial change of our findings that 
FM is the most stable phase. In any case, the structural degen- 
eracy for FM ZrMn2 would be unaffected. 

So far, the polymorphism of ZrMn2 is a predicted property 
based on calculated data. The question is how reliable are the 
results derived by the applied ab initio approach? An estab- 
lished fact is the occurence of the C14 structure at elevated 
temperatures, which can be explained by our theoretical ap- 
proach on taking into account the lattice vibrations. From the 
calculated phonon dispersion relations and phonon densities 
of states for C14 and C15 the temperature dependent free vi- 
brational energy Fyib{T) is derived. Based on the difference 
curves in Fig. |3lwe predict a structural phase transition CI 5 
^ C14 at a transition temperature of Ttr ~ 160K. A temper- 
ature dependent free energy contribution for the electrons was 
included as well fisll . which, however, is unimportant at T = 
Ttr. For temperatures larger than Tt^ C14 is stabilized at least 
in comparison to the CI 5 structure. (Because of its large unit 
cell we did not calculate the phonon dispersion relations for 
C36) Decomposing the free energy Fyn, = Umb — TSmb in 
contributions due to the internal energy Umb and the entropy 
Svib we find that the C15 ^ C14 transition is caused by the vi- 
brational entropy. Already the zero point energy Fyib (T = 0) 
contribution reduces the energy difference between C14 and 
C15. Although the C14 and C15 structures are nearly ener- 
getically degenerate no indication of any structural instability 
is seen in the phonon dispersion relations LI 41 as would be 
typical for martensitic transformations 0]. Finally, the elas- 
tic constants calculated for the CI 4 and C15 structures are 
typical for stable Laves phase compounds (Til and from the 
conditions for mechanical stability against any homogeneous 
elastic deformations Il9l such instabilities can be excluded for 
both structures. 
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FIG. 3: Difference of free energies F(C14) - _F(C15) of ZrMn2 in 
the ferromagnetic state. Shown are results with contributions of lat- 
tice vibrations only (full line) and also with including electronic con- 
tributions (dashed line). Positive values: C15 structure is preferred; 
negative values: C14 is stable. Insert magnifies region around the 
temperature Itr of the structural transition C15 C14. 



The reliablity of our calculational approach is strongly con- 
firmed by results on related compounds. For ZrCr2 12011 and 
ZrFe2 '21 J recent experiments find the C15 structure stable 
up to temperatures close to the melting point. When studying 
the three structures for these two compounds we also find that 
the C15 structure is definitely preferred |14]. Interestingly, 
according to our calculations ZrCr2 is nonmagnetic but ZrFe2 
is found to be strongly ferromagnetic. Analogous to ZrMn2 
the C15 structure gains significantly more magnetization en- 
ergy as compared to the hexagonal cases. Because for ZrFe2 
the magnetic moment of 1.9 jiB (derived from our calcula- 
tion as well as from experiment) is significantly larger than for 
ZrMn2 the magnetic energy gain of the C15 structure is also 
enhanced and C15 is now stabilized. For ZrMn2, as discussed 
above, the magnetization energy leads to the degeneracy of 
structures. For TiMn2 our ab initio approach yields the cor- 
rect C14 structure ground state 1 16]. The compound TiMn2 is 
just nonmagnetic (C14) or weakly magnetic (CI 5). 

An indicative quantity about the strength of bonding is the 
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2Fa{Mn) defined by the differences of free energies at the 
equilibrium volumes for the compound, Fcomp, and the pure 
elements, Fq, respectively. At zero temperature the expres- 
sion reduces to the differences of the corresponding energies 
Udft neglecting the weak temperature dependence for mod- 
erate temperatures sufficiently lower than the melting temper- 
ature. From the differences between the corresponding en- 
ergies, Udft, we arrive at a formation energy of about -64 
kJmol"^ for the ferromagnetic ground state, which character- 
izes a rather strong bonding between the Zr and Mn atoms. 
Therefore, from the formation energy no indication of any 
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phase instability can be inferred. The available experimental 
information of -48 kJmol^^ in Ref. |23] needs to be checked 
by more recent experiments. The strong formation energy of 
ZrMn2 is quite in contrast to the very closely related com- 
pound YMn2. From standard DFT calculations we derive 
even slightly positive formation energies, and experimentally 
the compound hardly forms | l4 l22ll . The reliablity of our cal- 
culated energies of formation is confirmed by the results for 
TiMn2 of about -90 kJmol^^, on which experiment and theory 
agree very well| 16]. 

For the C14 structure of ZrMn2, for which experimental 
data exists, the calculated structural and lattice parameters are 
in good agreement with experiment 1 14]. The calculated vol- 
ume is slightly smaller by about 3.4% possibly due to the ap- 
proximations made for the many-body term of DFT. It should, 
however, be noted that the deviations are small and that our 
basic findings of energy degeneracy involve differences of en- 
ergy, each of them obtained with the same approximations. 
Furthermore, the qualtity of our data is corroborated by the 
calculated energies of formation which in many cases agree 
with reliable experimental data. Very recently, ab initio DFT 
calculations for the C14 and C15 structures of ZrMn2 were 
published ]24]. We find perfect agreement for the magnetic 
moment but -due to the lack of published data- no comparison 
of structural energy differences is possible. 

Summarizing our results, based on density functional cal- 
culations we predict ZrMn2 to be polymorphic at low tem- 
peratures. We propose the degeneracy of at least three crys- 
tal structures (C14, C15, C36) which are distinguishable by 
different stacking of layers. This degeneracy occurs in the 
ferromagnetic phases, which are lowest in energy compared 
to nonmagnetic or selected antiferromagnetic spin orderings. 
The gain in magnetization energy due to spontaneous spin po- 
larization is different for all investigated structures but just of 
the right size to cause the degeneracy. The CI 5 structure is 
stable by a very small energy difference which might disap- 
pear when pressure is applied. Including the free energies of 
lattice vibrations we predict a structural phase transition at 160 
K to the C14 structure which is found experimentally at ele- 
vated temperatures. The structural degeneracy is not due to 
any instability: the phonon dispersions for the C14 and C15 
structures show no soft modes, the energy of formation indi- 
cates strong bonding, and the elastic constants for C14 and 
C15 exclude any weakness due to elastic deformations. Fi- 
nally, the reliability of our calculations is confirmed by com- 
parison to reliable experimental data, for related compounds 
as well. 
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